Dynamics in a supercooled liquid of symmetric dumbbells: 
Reorientational hopping for small molecular elongations 
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We present extensive molecular dynamics simulations of a liquid of symmetric dumbbells, for 
constant packing fraction, as a function of temperature and molecular elongation. For large elon- 
gations, translational and rotational degrees of freedom freeze at the same temperature. For small 
elongations only the even rotational degrees of freedom remain coupled to translational motions and 
arrest at a finite common temperature. The odd rotational degrees of freedom remain ergodic at all 
investigated temperature and the temperature dependence of the corresponding characteristic time 
is well described by an Arrhenius law. Finally, we discuss the evidence in favor of the presence of a 
type-A transition temperature for the odd rotational degrees of freedom, distinct from the type-B 
transition associated with the arrest of the translational and even rotational ones, as predicted by 
the mode-coupling theory for the glass transition. 
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The ideal mode-coupling theory (MCT) equations have 
been recently solved in the site-site representation for a 
system of symmetric hard dumbbell molecules Q, Q , as 
a function of the packing fraction <p and the elongation 
£. Interestingly enough, the theory predicts two different 
dynamic arrest scenarios, on varing tp and £ (see Fig. 1 
in Ref. pj). For large elongations, it is predicted that 
all rotational correlation functions are strongly coupled 
to the translational degrees of freedom and dynamic ar- 
rest takes place at a common tp value, ipf((^). According 
to MCT, the transition is of type-B, i.e. the long time 
limit of translational and rotational correlation functions 
jumps discontinuosly from zero to a finite value at the 
ideal glass transition line. The ideal glass transition line 
(the _B-line) has a non-monotonic shape in the tp— £ plane, 
with a maximum at £ = 0.43. The i?-line continues for 
small elongations until the hard sphere limit at £ = is 
reached. However, for small elongations ( < ( c = 0.345, 
theory predicts a novel different scenario: only the trans- 
lational and the even rotational degrees of freedom freeze 
at the B-line. Here even and odd refer to the parity of the 
order I of the rotational correlator Ci (t) = (Pj(e(t)-e(O))), 
where Pi is the Legendre /-polynomial and e(t) is a unity 
vector along the molecular axis at time t. The odd-l ro- 
tational degrees of freedom freeze at higher values of the 
packing fraction, namely at an A-line pf(() > <Pc(0> 
that merges with the £?-line at ( = £ c . The glass tran- 
sition of the odd degrees is of type-A, i.e. the long time 
limit of the odd rotational correlation functions increases 
continuously from zero on crossing the A-line. 

The A- and B-lines separate the tp — £ plane in three 
dynamic regions: an ergodic fluid, a completely arrested 
state, and an intermediate state (located between the A- 
and i?-lines) which is the amorphous analog of a plastic 
crystal. In such a state, each molecule remains trapped 
in the cage formed by the neighbouring molecules, but 



is able to perform 180° rotations within the cage, which 
lead to relaxation for the odd-l but not for the even-/ 
rotational correlators. An analogous scenario is predicted 
by molecular MCT for the case of hard-ellipsoids in the 
limit of small aspect ratio Q. 

In order to test these predictions, we have recently car- 
ried out molecular dynamics simulations in a binary 
mixture of Lennard- Jones (LJ) dumbbell molecules at 
fixed packing fraction, as a function of the elongation 
using the temperature T as control parameter (instead 
of tp used in the theoretical work). In agreement with 
MCT predictions, a non-monotonic dependence on the 
elongation has been observed for the isodiffusivity curves, 
with a minimum at a value of £ close to the maximum of 
tPciO M- For small elongations Q < 0.3, in deep super- 
cooled states, the coherent intermediate scattering func- 
tion evaluated at the maximum of the static structure 
factor, as well as the rotational correlator C^t), decays 
to zero at times several orders of magnitude longer than 
those for C\ (t) . An analysis of the plateau height for the 
C\ (t) correlator shows a sharp drop of the former in the 
range 0.3 < £ < 0.5, consistent with the theoretical £ c 
value predicted for the hard dumbbell fluid. These fea- 
tures suggest the presence of a nearby type- A transition 
in the investigated LJ dumbbell system. 

Simulations in Ref. 4] have been carried out in 
equilibrium liquid states, for state points above the B- 
transition, i.e., at T such that translational, and both the 
odd and even rotational degrees of freedom are ergodic. 
In this manuscript we complement equilibrium simula- 
tions and we also investigate the reorientational dynamics 
below the i3-transition, i.e., in out-of- equilibrium states, 
where the translational and even rotational degrees of 
freedom are arrested. We aim at providing evidence of 
the type-^4 transition, and to shed light on the dynamic 
features of amorphous analogs to plastic crystals. 
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FIG. 1: Isodiffusivity lines (thin continuous lines), and ti 
(dotted lines) and t 2 (thin dashed lines) isochrones in the T- 
C, plane. These quantities are averaged over all particles in 
the system. Line segments join the numerical data (symbols 
are not included for clarity). The corresponding values (from 
top to bottom) are: D = 1 x 10~ 2 ,3 x 10~ 3 , 1 x 10" 3 ,3 x 
10" 4 , 1 x 10" 4 , 3 x 1(T 5 ; ti = 3 x 10 1 , 1 x 10 2 , 3 x 10 2 , 1 x 
10 3 ,3 x 10 3 , 1 x 10 4 ,3 x 10 4 , 1 x 10 5 ; r 2 = 3 x 10 1 , 1 x 10 2 ,3 x 
10 2 , 1 x 10 3 ,3 x 10 3 , 1 x 10 4 . Empty and filled symbols are 
respectively the estimated MCT and VFT temperatures (see 
text) for D (circles), n (squares), and r 2 (triangles). Thick 
dashed and thick continuous lines are corresponding guides 
for the eyes. 

We investigate a binary mixture of 410:102 dumbbell 
molecules [B| in a box of side L with periodic bound- 
ary conditions. Each molecule consists of two identical 
atoms, denoted by A and B respectively for the most and 
the less abundant component of the mixture. The inter- 
action between atoms of different molecules is given by 
a Lennard- Jones potential plus a linear term: V a p{r) — 
^al3[(o a f3/r) 12 - (<J a p/r) 6 + A + Ai{r/a a p)], where a, 
j3 G (A, B). The interaction parameters are the same 
of the monoatomic binary mixture of Kob and Andersen 
@, i.e., e AA = 1.0, e AB = 1.5, e B B = 0.5, a AA = 1.0, 
cab = 0.8, obb = 0.88. The values of A and A\ are 
selected to guarantee continuity of potential and forces 
at the cutoff distance ca a p, with c = 2.5. The resulting 
values are A = c~ 6 (7 - 13c~ 6 ) and A l = 6c- 7 (2c~ 6 - 1). 
The same atomic mass m = 1 is used for both A and B 
particles (each dumbbell has mass equal to two). 

Simulations have been carried out for elongations ( = 
0.125, 0.20, 0.30, 0.40, 0.55, 0.70 and 0.85 at different 
values of T and constant packing fraction ip = 0.708. 
The packing fraction is defined by <p = (ir/6)(N AA a AA + 
N bb <j 3 bb )L- 3 [1 + (3/2)C - (1/2)C 3 ], where N AA = 410, 
Nbb = 102 are the number of AA and BB molecules. 
Temperature, time and diffusion coefficient D are re- 
spectively calculated in units of e AA /k B , (ma AA / 1 t AA ) 1 ^' 1 
and (a\ A e AA /m) 1 / 2 . Equations of motion are integrated, 



with a time step ranging from 2 x 10~ 4 to 5 x 10~ 3 ac- 
cording to T and/or (, using the velocity Verlet scheme 
combined with the SHAKE algorithm for keeping the 
molecular elongation constant. At each state point, the 
system is thermalized at the requested T by periodic ve- 
locity rescaling. A microcanonical run is performed after 
suitable equilibration to generate a representative trajec- 
tory in phase space. Statistical averages are performed 
over 20 independently generated equilibrium trajectories. 
At high T, the thermalization process is started from 
a distorted crystalline configuration that quickly melts. 
At low T, it is started from equilibrated high T con- 
figurations. After reaching equilibrium, energy or pres- 
sure show no drift. Also, translational and rotational 
mean squared displacements as well as dynamic correla- 
tion functions, show no aging, i.e. there is no systematic 
dependence of the dynamical properties from the selected 
time origin. In practice equilibration requires that each 
dumbbell, in average, has moved several diameters and 
performed several full rotations from its initial position. 

At very low T and small £, relaxation times for trans- 
lational an even rotational degrees of freedom become 
exceedingly long for the time scale of the simulation. In 
these conditions only the odd rotational degrees of free- 
dom decay to zero. We perform long runs of ~ 10 steps 
until neither significant drift in T and pressure, nor aging 
in the translational and even rotational correlation func- 
tions is observed, to ensure that the residual rotational 
motion of the dumbbells takes place in an essentially 
frozen structure 8J. After this preliminary simulation, 
trajectories are calculated. 

Fig. ^ shows data for D and for the relaxation times 
Ti , T2 of the rotational correlators C\ (i) and C-2 (t) . The 
diffusion coefficient is evaluated from the long time limit 
of the mean square displacement, while t\ and are de- 
fined as the time at which the corresponding correlation 
functions decay to 3 % of their initial value. As previ- 
ously done in the analysis of other systems characterized 
by slow dynamics |J, |9| , we report isochrones (constant 
Ti) and isodiffusivity (constant D) curves, to provide an 
indication of the shape of the glass line in the (T — () 
plane. The isodiffusivity curves are consistent with the 
same quantities calculated in Ref. Q and, as previously 
noted, show a mimimum around £ = 0.55, close to the 
MCT prediction for hard dumbbells 4J. Fig. ^ shows 
that, especially at low T, the iso-T2 curves follow the 
same trend as the isodiffusivity curves, tending to col- 
lapse in a common curve, as expected on approaching a 
type-B transition. The iso-ri curves show a similar trend 
only at large £■ O n the contrary, for small elongations 
the iso-Ti curves depart from the iso-D and iso-T2 curves. 
Differently from ti , within the resolution in £ of our cal- 
culations, the T\ isochrones are monotonic in the T — Q 
plane. The fact that the iso-Ti curves lay significantly 
below the iso-T2 curves shows very clearly that fast re- 
laxation of the odd degrees of freedom takes place even 
when the even degrees are arrested. 

A detailed quantitative comparison of the presented re- 
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FIG. 2: Symbols: Temperature dependence of D, ti and T2. 
Continuous lines are fits to VFT functions (see text). Dashed 
lines are Arrhenius fits. Data for n have been represented vs. 
l/T in order to evidence the small £ Arrhenius behavior. 

suits with MCT predictions for dynamic exponents and 
scaling laws can not be simply performed, except for large 
(. Indeed, it has been shown that the range of validity 
of the asymptotic MCT laws — for the small ( hard 
dumbbell fluid case Q — is limited to unusually small 
values of the separation parameter e = (ipf — ip)/if B . 
In particular, the exponent 7 for the well-known power 
law D, 1/t oc (ifc — ip) 7 can only be correctly estimated 
by exploring the range e < 10 . For small £, a fit 
of data in the range 10^ 2 < e < 1CP 1 would show an 
apparent power law behavior, but the corresponding fit 
would yield an "effective" exponent, which contrary to 
the asymptotic limit, would not take a universal value, 
but would depend on the chosen observable Q . A similar 
phenomenon is expected to be present in the studied sys- 
tem at small £, preventing the possibility of an accurate 
determination of the asymptotic laws close to the MCT 
critical temperature T C B . In addition, we recall that ideal 



MCT neglects activated processes which are known to be 
present in reality around and below T B | 10| . Only when 
hopping phenomena can be neglected, an MCT study can 
be feasible very close to T C B . 

Fig. [2] shows the T dependence of D, t\ and T2 for 
several £ values. If a fit (not shown) of the T dependence 
of D and I/V2 according to the functional form (T — T^ 3 ) 1 
is performed, we obtain very similar values of T c s for 
both quantities in all the £ range (see Fig. 0. For the 
case of 1 /t\ , reasonable power-law fits can be performed 
only for £ > 0.4. For ( — 0.55 fits provide values of 
T C B consistent with those obtained for D and I/V2 but 
with different 7 values. Only for £ > 0.70 a simultaneous 
consistent description for D, 1/t± and I/V2 is achieved, 
with identical T C B and 7 values. Specifically T C B = 0.732, 
7 = 2.50 for C = 0.70 and T C B = 1.02, 7 = 2.06 for 
C = 0.85. 

To provide a lower bound to T B at small £, we also 
fit the data in terms of the well-known phenomcno- 
logical Vogel-Fulcher-Tammann (VFT) law £>,1/ti2 oc 
exp[— A/(T — Tvft)], which predicts a dynamic arrest 
of the corresponding degree of freedom at a finite tem- 
perature T VFT < Tf (l(J. In all cases, the VFT law 
provides a consistent effective representation of the dy- 
namics, with a fit quality comparable to the MCT power 
law fit. The resulting Tvft(C) f° r D, t\ and T2 are re- 
ported in Fig. ^ Interestingly, Tvft for D, t\ and T2 
coincide for large elongations, indicating strong couplings 
between the translational and all the rotational degrees 
of freedom. However, for ( < 0.4 the Tvft for n is 
much smaller than the Tvft for D and t%, confirming 
that around this value of ( the odd degrees of freedom 
completely decouple form the even and from the transla- 
tional ones. This latter value of ( is consistent with the 
merging point of the A- and _B-transitions for the hard 
dumbbell fluid and with estimates of Ref. 0. At 
smaller elongations £ < 0.3, for n we find Tvft ~ 0. 
This implies that a good description of the data is pro- 
vided by an Arrhenius law n oc exp(i?/fcsT) (see Fig. [3J|. 
The activation energies E/caa = 1-8, 3.3 and 5.0, respec- 
tively for £= 0.125, 0.20 and 0.30, roughly follow a linear 
dependence on £. The presence of an activated Arrhenius 
T dependence in T\(T) for small molecular elongations, 
strongly suggests that the residual reorientational dy- 
namics at very low temperatures, in the T-region where 
the translational dynamics is already arrested, is charac- 
terized by 180° rotations which do not lead to relaxation 
of the even rotational degrees of freedom. The Arrhenius 
functional form also suggests that hopping phenomena 
are extremely relevant, so much to completely mask the 
presence of a finite type-^4 glass transition temperature, 
if the analysis is limited only to the T dependence of n 
(or any other odd r). 

While we can not expect, due to this intense hopping 
decorrelation, to observe any critical law, indication of 
a residual cross-over finite temperature can still be ob- 
served in the full dependence of the correlation functions, 
and in particular in the strength of the long time a relax- 
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FIG. 3: Temperature dependence of the rotational correlators 
C<i(t) and C3(t) vs. time rescaled by T 1 ^ 2 , for elongation 
C = 0.2. From top to bottom T= 0.217, 0.276, 0.343, 0.496, 
0.750, 1.110, 1.396, 1.693, 1.770, 1.815, 1.933, 2.443 and 3.776. 

ation process. The main effect of hopping is indeed man- 
ifested in the long time decay from the plateau value. In- 
deed, in the absence of these additional relaxation chan- 
nels, the correlation function would decay only down to 
the plateau value. Fig. [21 shows, for £ = 0.2, the In- 
dependence of C*2(i) and C^{t), in a wide range of T, 
(across Tvtf = 0.96 for T2). Though C\(t) has more 
experimental relevance, C 3 (t) is more adequate to dis- 
cuss the behavior of the non-ergodicity parameter, since 
it shows a more pronounced plateau. A representation 
of Ci(t) vs. tT 1 / 2 , i.e rescaling the time by the thermal 
velocity, yields a clearer visualization of the decay of the 
correlation function. Data reported in Fig. [3] show two 
important features: 

(i) at low T the plateau height (i.e., the non-ergodicity 
factor) for both Ci and C3 progressively grows up, signal- 
ing a crossing of the critical MCT temperatures. How- 
ever, while in case of C2(t), the increase of the plateau 



value is connected to a non-ergodic long time behavior 
(i.e. C2 does not decay to zero within the simulation time 
window), Cz{t) always relaxes toward zero. This fact 
evidences that hopping processes restore ergodicity for 
reorientations at T far below the temperature at which 
the plateau value starts to increase from zero, (ii) The 
decay of the C2 correlation function is characterized by 
a finite T interval in which the plateau value is essen- 
tially constant (within the present numerical error). For 
smaller T, the height of the plateau starts to increase, as 
discussed above. Instead, the decay of the C3 correlation 
function is not characterized by any plateau at high T, 
and the increase of the plateau starts from zero, consis- 
tently with the presence of a type- A transition which, due 
to the hopping effects, is significantly weakened. More- 
over, the first T at which the C3 plateau value is different 
from zero is lower than the T at which the C2 plateau 
starts to increase. This suggests that the type- A transi- 
tion temperature is different and lower than the type-i? 
T C B . We also note that the a-relaxation process for odd 
Ci is highly stretched, which might reflect the local cage 
heterogeneity. 

In summary, we have provided evidence for the exis- 
tence of the type-A transition predicted by MCT for the 
reorientational dynamics of glass-forming liquids of small 
molecular elongation. Our numerical results complement 
those recently reported for T > Tf in Ref. 4]. By 
computing relaxation times and plateau heights of odd 
rotational dynamic correlators — in the arrested state 
for the even rotational and the translational degrees of 
freedom — we show that reorientational hopping domi- 
nates the dynamics and restores ergodicity down to tem- 
peratures far below the type-A transition. Though these 
results have been specifically obtained for a liquid of sym- 
metric dumbbells, experimental techniques probing C\ (t) 
and Ciit), at very low temperatures, will yield qualita- 
tively similar results for small elongated molecules with 
a weak asymmetry which disfavours rotations different 
from ss 180°. Finally, the present model can become a 
test case for developing theories of the glass transition 
in which activated processes — which become relevant 
around and below — can be studied. 
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